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Let B — {Bi{t), . . . , Bdit)) be a d-dimensional fractional Brownian motion 
with Hurst index a < 1/4, or more generally a Gaussian process whose paths 
have the same local regularity. Defining properly iterated integrals of B is 
a difficult task because of the low Holder regularity index of its paths. Yet 
rough path theory shows it is the key to the construction of a stochastic 
calculus with respect to B, or to solving differential equations driven by B. 

We intend to show in a forthcoming series of papers how to desingularize 
iterated integrals by a weak singular non-Gaussian perturbation of the Gaus- 
sian measure defined by a limit in law procedure. Convergence is proved by 
using "standard" tools of constructive field theory, in particular cluster ex- 
pansions and renormalization. These powerful tools allow optimal estimates 
of the moments and call for an extension of the Gaussian tools such as for 
instance the Malliavin calculus. 

This first paper aims to be both a presentation of the basics of rough path 
theory to physicists, and of perturbative field theory to probabilists; it is only 
heuristic, in particular because the desingularization of iterated integrals is 
really a non-perturbative effect. It is also meant to be a general motivating 
introduction to the subject, with some insights into quantum field theory 
and stochastic calculus. The interested reader should read in a second time 
the companion article |48] (or a preliminary version [47]) for the constructive 
proofs. 
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Introduction 



A major achievement of the probabilistic school since the middle of the 20th 
century is the study of diffusion equations, in connection with Brownian mo- 
tion or more generally Markov processes - and also with partial differential 
equations, through the Feynman-Kac formula - with many applications in 
physics and chemistry [64] . One of the main tools is stochastic calculus with 
respect to semi-martingales M. An adapted integral such as X{u)dM{u) 
may be understood as a limit in some sense to be defined. Classically one 
uses piecewise linear interpolations, Y2s<ti<...<tj^<t 

X{ti)[M{ti+i) - M{ti)) 

or 'Yl,s<ti< <tN<t "^^*'-'"'"^^*'+^^ (M(tt-i-i) — M{ti)); these approximations de- 
fine in the limit N ^ oo the Ito, resp. Stratonovich integral. The latter one 
is actually obtained e.g. if M = W is Brownian motion and X{t) = f{Wt) 
with / smooth as the limit Hui^^q f{Ws{u))dWi;{u) for any smooth ap- 
proximation {We)e>o of W converging a.s. to W (see [69], or [35] p. 169). 
The Stratonovich integral X {u)d^^^"'^° M (u) has an advantage over the Ito 
integral in that it agrees with the fundamental theorem of calculus, namely, 
F{M{t)) = F{M{s)) + jlF'{M{u))d'^*'-''*°M{u). 

The semi-martingale approach fails altogether when considering stochas- 
tic processes with lower regularity. Brownian motion, and more generally 
semi-martingales (up to time reparametrization) , are (l/2)~-Holder, i.e. a- 
Holder for any a < 1/2 0. Processes with a-Holder paths, where a <C 1/2, 

^Recall that a continuous path X : [0, T] — >■ R is a-Holder, a £ (0,1), if 
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are maybe less common in nature but still deserve interest. Among these, 
the family of multifractional Gaussian processes is perhaps the most widely 
studied [52], but one may also cite diffusions on fractals [32], sub- or su- 
perdiffusions in porous media [261 [37| and the fascinating multi-fractal ran- 
dom measures/walks in connection with turbulence and two-dimensional 
Liouville quantum gravity [9l [13] . Many models in hydrodynamics take as 
input a space-time noise which is often chosen colored in space [34] . In this 
respect, let us mention in particular the Kraichnan model for passive advec- 
tion of scalars, for which anomalous correlation exponents [36[ [T^ [6] may 
be expanded in a for a — t- 0. 

We concentrate in this article on multiscale Gaussian processes (the ter- 
minology is ours) with scaling dimension or more or less equivalently Holder 
regularity a G (0,1/2), the best-known example of which being fractional 
Brownian motion (fBm for short) with Hurst index a, B'^{t) or simply B{t) 
H. We consider more precisely a two-dimensional fBm, B{t) = {Bi{t), B2{t)), 
with independent, identically distributed components El- The covariance ker- 
nel 'KBi{s)Bj{t) = ^5jj(|sp"-|- — |t — sp°) is that of an integrated colored 
noise in the physical terminology 0. It is a process with long-range, negative 
correlations, which is quite unusual from a statistical physics point of view; 
but the emphasis here is on the short-distance (or ultra-violet) behaviour, 
not on the long-distance one. 

The simplest non-trivial stochastic integral is then 

A{s,t):= I dBiiti) rdB2{t2)= I {B2{u) - B2{s))dBi{u), (0.1) 

a twice iterated integral, where B = {Bi{t), B2{t)) is a two-component fBm 
with independent, identically distributed components. Since 

/ dBi{ti) / dB2{t2)+ / dB2{t2) / dB^iti) = {Bi{t)-Bi{s)){B2{t)-B2{s)), 

J S 'J S o S i/S 

one is mainly interested in the antisymmetrized quantity (measuring a signed 

^It is (up to a constant) the unique self-similar Gaussian process with stationary incre- 
ments. The last property implies that its derivative is a (distribution-valued) stationary 
field. 

^The one-dimensional case is very different and much simpler, and has been treated in 

m- 

*In other words (informally at least) 'E.B[{s)B'j{t) ~ —Ca\t — sp" instead of 5{t — s). 
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follows from the Green-Riemann formula), 



£A{s,t) := I dBi{ti) r dB2{t2)- j dB2{t2) r dBi{ti) 

t 

{B2{u) - B2{s))dBi{u) - {Bi{u) - Bi{s))dB2{u), 

(0.2) 

called Levy area. The corresponding Stratonovich mtegral, obtained as a 
limit either by linear interpolation or by more refined Gaussian approxima- 
tions im [5H [58l 159] . has been shown to diverge as soon as a < 1/4. 

This seemingly no-go theorem, although clear and derived by straight- 
forward computations that we reproduce in short in section 1, appears to 
be a puzzle when put in front of the results of rough path theory |431 HH 
[28t [39} [20] . The essential idea conveyed by this theory - we shall make 
this precise in section 2 - is that a path F : M — )• with Holder regularity 
index a G (0, 1) must be seen as the projection onto the d first components 
of some "essentially arbitrary" rough path over T, denoted by 

r : M2 9 (s,t) ^ Tts := {Vl, . . . , T^) G M'^ x M"^" x . . . x M'^^, (0.3) 

which may be interpreted as iterated integrals of F in a limiting 
sense, namely, lim^^^o drf (ti) = Fj^ (t) — F,^ (s) for the d first components 
and 

/t f-ti pt rtN-i 

dTUti)j^ dTUt2),...,limJ^dTlit^)...J^ dTlitn) (0.4) 

for the remaining ones, for some smooth family of approximations (F^)£>o 
of F. The limit must be understood in a Holder norm sense, as explained 
in section 2. In other words, there exist infinitely many different families 
of approximations of B leading to as many different definitions of its iter- 
ated integrals! Alas, Gaussian approximations are unfortunately seemingly 
unable to produce such a definition for fBm with Hurst index a < 1/4. 

Our project in this series of papers is to define an explicit rough path 
over fBm with arbitrary Hurst index, or more generally multiscale Gaussian 
fields (see a companion article |48j and [Uj for a preliminary version ) by 
an explicit, probabilistically meaningful limiting procedure, thus solving at 



^where [ . J stands for the integer part of its argument. 
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last the problem of constructing a full-fledged, Stratonovich-like integration 
with respect to fBm. 

Let us explain our strategy for 1/6 < a < 1/4. Roughly speaking, our 
rough path is obtained by making B = (-B(l), .8(2)) interact through a weak 
but singular quartic, non-local interaction, which plays the role of a squared 
kinetic momentum, or bending energy, and makes its Levy area - and at 
the same time the iterated integrals of higher order - finite. Following the 
common use of quantum field theory, this is implemented by multiplying 
(probabilists would say: penalizing) the Gaussian measure by the exponen- 
tial weight e-3^'-i'/^-*(<^i'<^2)(ti,t2)|ti-t2|-^"dtidt2 [e] ^ ^^^^i 

Ani(0i,02)(ti,t2) = {{dA+){h){dA+){t2) + idA-){h){dA-){t2)} , 

(0.5) 

where: A (the coupling parameter) is a small, positive constant; (j)i,(j)2 are 
the (infra-red divergent) stationary fields associated to i?i,i?2, with covari- 
ance kernel as in eq. (II. 4p . and similarly, A^ are stationary left- and right- 
turning fields, built out of cpi , 02 and representing the singular part of the 
Levy area (see section 1 for details). As usual in quantum field theory, one 
considers first the truncated measure obtained by an "ultra-violet cut-off^' 
and on a finite "volume" (or finite horizon, in the probabilistic terminology) 
V = [—T, T], i.e. one multiplies the Fourier transforms of the fields 4>i,4>2 by 
some cut-ofi^ function with compact support in [—M^, ilf] (for some fixed 
constant M > 1) and integrates over V; see Definition 13.11 for the precise 
procedure. Then dA^ are replaced by the truncated quantities {dA^)^^ 
built out of the truncated fields The truncated interacting Lagrangian 
reads 



ic'^ [ [ \tl-t2\-^''CZ?{(pl,(t>2){tl,t2)dtidt2+ [ C 
^ J JvxV Jv 



P 

bdry 



kA'/ / \h-t2\-"'{{^A+rP{t^)i^A-^rPit2) 

^ J JvxV 



+{dA-rf{t,){dA-)^p{t2)}dtidt2 + J C 



p 

bdry' 



(0.6) 



where is some singular "Fourier boundary term" multiplied by an 

evanescent factor M~'^^ (k > 0), which cures unwanted difficulties due to 



®The unessential constant c'^ is fixed e.g. by demanding that the Fourier transform of 
the kernel c^jfi - 42]"*° is the function 



5 



the ultra-violet cut-off 0. When p and V are finite, the underlying Gaussian 
fields are smooth, which ensures the existence of the penalized measure. The 
assertion is that the penalized measures converge weakly when p,\V\ — > oo 
to some well-defined, unique measure, while the truncated iterated integrals 
themselves converge in law to a rough path over B. 

Note that the statistical weight is maximal when dA'^ = dA^ = 0, i.e. 
for sample paths which are "essentially" straight lines. Another way to mo- 
tivate this interaction (following an image due to A. Lejay) is to understand 
that the divergence of the Levy area is due to the accumulation in a small 
region of space of small loops [39]; the statistical weight is unfavorable to 
such an accumulation. On the other hand, the law of the quantities in the 
first-order Gaussian chaos, characterized by the n-point functions 



{Bi^{xi)...Bi„{Xn))x 



ii, . . . , i„ = 1, 2, where 
Z ■.= E 



e 2 



c'q / f Ci„ti<l>l,(p2)itl,t2)\tl-t2\-'^"dtldt2 



(0.7) 



(0.8) 



is a normalization constant playing the role of a partition function, is in- 
sensitive to the interaction Thus we have built a rough path over fBm. 
This conveys the idea that the paths have been straightened by removing in 
average small bubbles of scale M~p. In doing so, the paths of the limiting 
process when p — ?> oo are indistinguishable from those of B, hut higher-order 
integrals have been corrected so as to become finite. 

Starting from the above field-theoretic description, the proof of finiteness 
and Holder regularity of the Levy area for A > small enough follows, 
despite some specific features, the broad scheme of constructive field theory, 
see e.g. the monographies [T| 149^156]. Constructive field theory is a program 
originally advocated in the sixties by A. S. Wightman [66], the aim of which 
was to give explicit examples of field theories with a non-trivial interaction 
; see Glimm and Jaffe's book [23j for an introduction and references therein 



^The exact form of C^J^y requires detailed constructive explanations and will not be 
required here. It is to be found in the companion article [48]. 

*In the two preceding equations, E [ ■ g-^'^a //^'"t('*i''*2)(ti.t2)|ti-t2r*°dtidt2 j stands 

for the limit of E [ ■ e'^'^'c //vxv £^?(^i.^2)(ti,t2)|ti-t2r"°dtidt2+/v j „hen ^ 
oo as we explained above. 
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for an extensive bibliography. Let us give a short guide to the history of the 
subject. 

The first contribution was made in 1965 by E. Nelson who introduced a 
scale analysis [50 1 to control the divergence of a model whose only divergence 
comes from Wick ordering. J. Glimm and A. Jaffe introduced the phase 
space analysis [23] for models having a finite number of divergent graphs. 
The cluster expansion was devised by J. Glimm, A. Jaffe and T. Spencer 
[25 j to control infinite volume limits. 

The Roman team [5] realized that the above phase space analysis was 
in some sense a continuous space version of the block-spin expansion, first 
written by Kadanoff for the Ising model, and then made into a major tool 
both in high-energy and statistical physics by K. G. Wilson through the 
introduction of the concept of renormalization group [67l |68] . The multi- 
scale expansion was devised in the eighties to provide a rigorous version of 
Wilson's renormalization group e.g. including the flow of the effective pa- 
rameters: see |22j for the block-spin approach, and [T7| for the continuous 
space multi-scale cluster expansion. 

For some fermionic theories a simpler version of these constructions is 
available, due to the fact that (contrary to the bosonic case) the series 
expansion in terms of the effective coupling constants is convergent [42 1 1^ . 
The multi-scale cluster expansion has also allowed to study models with a 
singularity around a surface, like the so-called jellium model of interacting, 
non-relativistic fermions \18\ I12j. modelling the generation of Cooper pairs, 
in connection with the famous BCS (Bardeen-Cooper-Schrieffer) theory of 
supraconductivity p9j . 

In this work we use the multi-scale cluster expansion developed in [17] 
more than twenty years ago which seems to us the most appropriate for 
these probabilistic models; it reduces to the minimum the use of abstract 
combinatorial identities and algebra, to the benefit of a very intuitive and 
visual (though sometimes heavy) tree expansion. 

The main theorem may be stated as follows. As a rule, we denote in this 
article by E[...] the Gaussian expectation and by (...)a,v,p the expectation 
with respect to the A-weighted interaction measure with scale p ultraviolet 
cut-off restricted to a compact interval V, so that in particular E[...] = 

(•••)o,oo- 

Theorem 0.1 Assume a S (^j^)- Consider for A > small enough the 
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family of probability measures (also called: {(p,d(l),a)-model) 
Fxy,p{(piA2) = e-^'''^^ ^'^'''^'^^''-'^^^'"^^-'^^''^^'''^''''^^^ 

(0.9) 

where dfi'^P{(j)i) = d^{(f)^ ^) is a Gaussian measure obtained by an ultra- 
violet cut-off at Fourier momentum [^| ~ (M > 1), see Definition \3.1[ 
Then (Pa, y,p)v,p converges in law when \V\,p — >■ oo to some measure ¥x, 
and the associated iterated integrals 

ft ftl ft ftl ftn-l 

d4>r/{ti) d<t>-:'it2),..., d^r/ih) d<p^/{t,)... d<pz'{tn),... 

Js J s J s Js J s 

converge in law to a rough path over B. 

The result is not difficult to understand heuristically, at least for quan- 
tum field theory experts, if one resorts to the non-rigorous perturbation 
theory (see sections 3 and 4). First, by a Hubbard-Stratonovich transforma- 
tion (a functional Fourier transform), one replaces the non-local interaction 
^{(f'lT 4'2){ti,t2)\ti — t2\~^°' with a local interaction C{4>i,(j)2,(T)(t) depend- 
ing on a two-component exchange particle field a = (cr+ (t), iT_(i)). Then 
a Schwinger-Dyson identity (a functional integration by parts) relates the 
moments of A to those of a. Simple power-counting arguments show that a 
connected 2n-point function of a alone is superficially divergent if and only 
if 1 — Ana > 0. Thus, restricting to a > 1/8, one only needs to renormalize 
the two-point function. Since the renormalized propagator of a is screened 
by a positive, infinite mass term, the theory is free once one has integrated 
out the (T-field, hence one retrieves the underlying Gaussian theory {(pi, 02). 
The Schwinger-Dyson identity then shows that the two-point functions of A 
have been made finite. Finally, simple arguments (not developed here) yield 
the convergence of higher-order iterated integrals in the interacting theory 
provided a > 1/6. 

Whereas these heuristic arguments are not difficult to follow in principle, 
they do not constitute at all a proof. Theorem 0.1 is proved in the com- 
panion article |48] by following - as explained above - the general scheme 
of constructive field theory. Although the constructive method is really a 
multi-scale refinement of the previous arguments, explaining it precisely is 
actually a formidable task, which is in general very much model-dependent, 
whereas perturbative renormalization always follows more or less the same 
lines; briefiy said, the difference lies in the difference between a formal power 
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series expansion and an analytic proof of convergence for a given quantity. 
This task we perform at long length and in great generality in the compan- 
ion article, with the view of making constructive arguments into classical 
mathematical tools which probabilists may eventually reemploy. 

Here is an outline of the article. 

We begin in Section 1 by recalling classical arguments (due to the second 
author) explaining the divergence of the Levy area for a < 1/4, which is the 
starting point for all the story [58]; Fourier normal ordering |62i [6T] - an 
indispensable tool for the sequel - is introduced there. Section 2 is a brief 
introduction into rough path theory, mainly for non-experts. Subsections 2.1 
and 2.2 are standard and may be skipped by experts, whereas subsection 
2.3 - a brief summary of the previous contributions of the second author to 
the subject - gives the context in which this series of papers arose. 

The heart of the article is Section 3 and Section 4. Our problem is recast 
into a quantum field theoretic language in section 3; we take the opportunity 
to explain the basis of quantum field theory and renormalization at the same 
time. The interaction term is introduced at this point, where it comes out 
naturally. Finally, section 4 is dedicated to a heuristic perturbative "proof 
of the convergence of the Levy area of the interacting process, and serves 
also in some sense as an introduction to the companion paper |48j . 



1 A Fourier analysis of the Levy area 

The quantity we want to define in the case of fractional Brownian motion is 
the following. 

Definition 1.1 (Levy area) The Levy area of a two-dimensional path T : 
M — 7- between s and t is the area between the straight line connecting 
(ri(s),r2(s)) to {Ti{t),T2{t)) and the curve {{Ti{u),T2{u)); s <u<t}. It 
is given by the following antisymmetric quantity, 

dv^iti) / dr2{t2) - / dr2{t2) / dViiti). (1.1) 

The purpose of this section is to show by using Fourier analysis why the 
Levy area of fBm diverges when a < 1/4. This is hopefully understandable 
to physicists, and also profitable to probabilists who are aware of other proofs 
of this fact, originally proved in [11], because Fourier analysis is essential 
in the analysis of Feynman graphs which shall be needed in section 4. We 
follow here the computations made in |61j or [60]. 
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Definition 1.2 (Harmonizable representation of fBm) LetW{(,),^ € 
be a complex Brownian motion H such that W{-0 = -W{C), and 

Bt := i27rc^y-2 teR. (1.2) 

J-oo K 

The field t E M is called fractional Brownian motion^^ . Its paths are 
almost surely a~ Holder, i.e. (a— e)-H61der for every e > 0. It has dependent 
but identically distributed (or in other words, stationary) increments Bf—Bs- 
In order to gain translation invariance, we shall rather use the closely related 
stationary process 

'^W ■= / -r\^\"''dWiO, tEM (1.3) 

- with covariance 

- which is infrared divergent, i.e. divergent around ^ = 0. However, the 
increments (j)(t) — 4'{s) = Bt — Bs are well-defined for any (s,t) E M^. 

In order to understand the analytic properties of the Levy area of fBm, 
we shall resort to a Fourier transform. One obtains, using the harmonizable 
representation of fBm, 

A{s,t) := J dBi{ti) dB2{t2) 

The Levy area CA[s, t) := CAb{s, t) is obtained from this twice iterated 
integral by antisymmetrization. Note that CA{s, t) is homogeneous of degree 
2a in \t — s\ since B[ct) — B{cs), c > has same law as c°'{B{t) — B{s)) by 
self-similarity. 

Expanding the right-hand side yields an expression which is not homo- 
geneous in ^. Hence it is preferable to define instead the following stationary 

^Formally, {W'i^i)W'i^2)} = and {W'{^i)W^) = 5{^i - ^2) if 6,^2 > 0. 
^"^The constant Ca is conventionally chosen so that E(_Bt — Bs)^ = 1* ^ sp". 
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quantity called skeleton integral, which depends only on one variable, 



A{t) := J dBi{ti) dB2{t2) 

1 [ dWi{ii)dW2{i2) 



t Hi 

dti / dt2 ■ e^(*i«i+*2«2) 

(1.6) 



2vrc„7 |eih-V2|^2|a-i/2 
1 f dWi{^i)dW2{^2) e'*(5i+«2) 



27TCaJ [i(ei+6)][ie2]' 

where by definition J^e'^^dn = From A{t) and the one-dimensional 
skeleton integral 

Mt) = (2vrc,)-5 J dB,{u) = j . —, (1.7) 

which is the above-defined infra-red divergent stationary process associated 
to B, one easily retrieves A{s, t) since 

A{s,t) = J^Biiti) dB2{t2) - J' dB2{t2)^ 

= A{t) - Ais) + Adis,t), (1.8) 

where {2TTCa)^ Ad{s,t) := (-Bi(t) — Bi{s))(p2is) (called boundary term) is a 
product of first- order integrals. 

One may easily estimate these quantities in each sector ^ |^2|- In 
practice, it turns out that estimates are easiest to get after a permutation 
of the integrals (applying Pubini's theorem) such that (for twice or multi- 
ple iterated integrals equally well) innermost (or rightmost) integrals bear 
highest Fourier frequencies; this is the essence of Fourier normal ordering 
[62 \ 116 ^ 163]. This gives a somewhat different decomposition with respect to 
(fLHIl since dBi{ti) f^^ dB2{t2) is rewritten as - Jl dB2{t2) dBi{ti) in 
the "negative" sector |^i| > |^2|- After some elementary computations, one 
gets the following. 

Lemma 1.3 Let 

^+(t) := 27TC^l'dti y*'dt2-F-l((6,6)^l|6|<|6|(-^^l)(6)(^52)(6))(tl,t2) 



dWi{^i)dW2ie. 



2 



[i(6 + 6)][i6] 



(1.9) 
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and 



A-{t) := 27rcal dt2y"'dii^-H(6,6)^l|6|<|«i|(-^^i)(6)(-^52)(6))(ii,i2) 
4l<l6i ■ [i(6+6)][i6]' ^' ^ 



Then 

As,t) = ^ {iA+{t) - A'is)) - {A-{t) - A-{s)) + {At{s,t) - As{s,t))} , 

(1.11) 

the boundary term Aq — Aq being given by 

At{s.t)-A-s{s,t) = 1- ^ 

(git6 _ eis6)ei*€i 1 dVFi (6)^^^/2(6) 



(1.12) 



Two lines of computations show immediately that 
Var^±(s,t) < ||ei*«-e-«|2|e|-i-^"dC 



< (1.13) 

so that (essentially by the Kolmogorov-Centsov lemma, see section 2) the 
Holder regularity indices of Bi and B2 add in the case of the boundary term, 
to produce a quantity which is 2a~-Holder. (Note that the artificial infrared 
divergence at 6 = disappears when Taylor expanding e'*^i — e'^^^). On 
the other hand, letting ^ := 6 + ^2 and introducing an ultra-violet cut-off at 
161 = A 1, one may see for instance A'^{t) as an inverse random Fourier 

transform of the integral ^ ^ i|g-6l<l6l '^'^if^^ |g_g^|a-i/2|g^|a-i/2 , whose 

variance diverges like ^ = 0(A^~^") or O(lnA) in the ultra-violet limit 

^2 

A ^ 00 as soon as a < 1/4. Note that the ultraviolet divergence is in the 
region |6|,|6| > 1^1 • 

It is apparent that the central role in this decomposition is played by the 
Fourier projection operator Z)(l|^j|<|^2|) = -^"^ (-'-|?i|<|42| ' *^( • )) • Since 
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are obtained by Fourier projecting {Bi{t) — i?i(s))02(s), or (i?2(i) — 
S2(s))(/>i(t), which are perfectly well-defined products of continuous fields 
PI . it was clear from the onset that these would be regular terms. Hence 
singularities come only from the one-time quantity A^{t), which does not 
split into a product of first-order integrals, and that we shall call the singular 
part of the Levy area. 



2 An introduction to rough paths 

2.1 General issues 

Let r = (ri(t), . . . ,Tfi{t)) be a smooth path with d components. As ex- 
plained in the Introduction, the Levy area of F, Tfg{i, j) := dTi{ti) f^^ dTj{t2) 
is the simplest non-trivial iterated integral of T. The interest for iterated 
integrals of T comes from the study of two closely related problems in the 
case when T is not regular any more. 

1. Integration along an irregular path. 

Assume one wants to define the integral of the (say, smooth) one-form 
/ := Yl'j=i fj{x)dx^ along the path F, namely, the quantity 

fdV := J2j Is fj(^i''^))d^j{'^)- Since F is not differentiable, drj{u) may 
not be understood as • du, and the very meaning of this quantity is un- 
clear. Unfortunately, the Riemann-type sum Yl7=o fj(^i'^i))(^ji'^i+i) ~ 
Fj(ij)), with s = to < . . . < ti = s + ^(t — s) < . . . < tn = t, may be shown 
to diverge in general as soon as a < ^ o. 

A Taylor expansion to order of the integrand yields (coming back to 
the case of a regular path) the improved Riemann-type sum 

n-l N d ^p_i „ 

EE E 9,. t- (r(t.))rf^,,, (,„..., J,), (2.1) 

i=0 p=lii,...jp=l ^i'" 



where 



K+.,uih^ ■■■Jp)= [ dTj^it^) ... f dT,^{tp) (2.2) 

J s J s 



apart from the spurious infra-red divergence (see above) 
^^From a naive bound by the 1-variation of the path, X]j X^i l^j (ij+i) — ~ 
0(n^~"), one would come to the erroneous conclusion that the Riemann-type sums di- 
verge when Q < 1. The so-called Young theory of integration (see e.g. [S^) lowers the 
barrier to a = i by taking into account the Holder regularity of the integrand f{T{t)). 
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is a p-th order iterated integral. The problem is, if F is irregular, iterated 
integrals of T are a priori ill-defined for the same reasons as before. 

2. Solutions of differential equations driven along an irregular 
path 

Consider the differential equation 

d 

dyt = Y,VMt))drj{t). (2.3) 

i=l 

The following series gives a formal solution, 

oo 

yt = ys+J2 E [Vn---Vi^-ld]{Ys) ■ T''{h,...,iN), (2.4) 

N=l l<ji,...,jjv<d 

with r as in eq. (|2.2p . Solutions are usually computed by using some 
iterated numerical scheme. For instance, the Euler scheme of rank N gives 
the solution to (j2.3p as the limit when n ^ oo of the compound mapping, 

HXt,t„_, ; • • • H^t,,t^ ; H^'"'; Vs)---), (2.5) 

where ^{Tts',ys) is the series (|2.4p truncated to order N. If one takes for 
F an a- Holder path, one stumbles again into the same problem of defining 

Tts = (r^s, . . . , Fj^). 

In both cases, the hope is that, if one finds some (non necessarily unique!) 
way of defining iterated integrals of F with the correct regularity properties, 
then the refined Riemann-type sums (j2.ip or Euler scheme (j2.5p converge 
when the mesh goes to 0. Rough path theory shows this is possible |^ 
provided one chooses > [1/aJ - here we choose = [1/aJ minimal - 
and 

Tts = {Ti{ii)i<i,<d, . . .,T^sih, . . .,^N)l<n,...,^M<d) ■ ^ M'^xM'^' x . . . xM' 

(2.6) 

is a rough path with Holder regularity index a in the following sense: 

Definition 2.1 (rough path) An a- Holder continuous rough path over F 
is a functional Ffg{ii, . . . , i„), n < [N := 1/qJ , ii, . . . , i„ G {1, . . . , d}, such 
that Tts{i) = Tt{i) — Ts{i) are the increments of T, and the following 3 
properties are satisfied: 

^^Furthermore, the hmit is a-H61der and satisfies nice continuity properties with respect 
to the path T. 
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(i) (Holder continuity) r"^(ii, . . . ,2„,) is na-Holder continuous as a func- 
tion of two variables, namely, sup^^jgjg * |t_g|a' ' < oo. 

(a) (Chen property ) 

(ii, ...,in)= r^uiil, ...,in)+ K,{ii, ...,in) + 

r'l^{il,...,in,)Kl(in,+l,...,in); (2.7) 

ni+n2=n 

(Hi) (shuffle property) 

rfs^C^l, • ■ ■ ,ini)rts^(jl, • • • ,in2) = rj^"^"*" ^(/ci, . . . , /Cn^+nj), 

k&Sh{i,j) 

(2.8) 

where Sh{i,j) - the set of shuffles of the words i and j - is the sub- 
set of permutations of the union of the lists i,j leaving unchanged 
the order of the sublists i and j. For instance, Tfg{ii,i2)Tlg{ji) = 
Tl{ii,i2,ji) + Tl{ii,ji,i2) + Tl{ji,ii,i2)- 

A formal rough path over T is a functional satisfying all the above prop- 
erties except Holder continuity (i). 

In a random setting, the Holder continuity estimates (i) are generally 
proved as a consequence of moment estimates such as E|r"^pP < Cp\t — 
g]2pna^ p > l^fi = _ _ ^ ]S[_ This may be seen as a consequence of the well- 
known Kolmogorov-Centsov lemma stating that (for a measurable process 
random T) 



{E[\T{t) - r(s)|2p] < C\t - ^ Va' < a, E 



irW-rMiy- 

:..€[0T| It-S\'- ) 



< oo 



(2.9) 

or more precisely of an extension (or a variant) of these estimates adaptated 
to functions of two variables (such as (s, t) i— t- T"g) due to Garsia, Rodemich 
and Rumsey [21]. 

In particular, if T is smooth, then its natural iterated integrals 
dTi^{ti) . . . dTi^{tn) satisfy properties (ii) and (iii). 

However, it is not clear a priori in what sense abstract data as in Defi- 
nition [2?T] should represent iterated integrals in the usual sense. 
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2.2 Geometric approach 

The answer to this question comes from a reinterpretation of rough paths 
in terms of group theory and geometric structures. We generally refer to 
the book by P. Priz and N. Victoir [2U] for this paragraph. Consider the 
signature Tts = (I^ts ' ^ts > • • • ) of ^ smooth path F as 

Tis, t):=l + Y^ Tiih)X^^ + Tlih,i2)X^^ CS> X^^ + . . . , (2.10) 

sitting inside the tensor algebra TM*^ = ©„>or"'M'^, with X^, . . . ,X'^ gener- 
ating a basis of M'^ ~ T^R'^. Note that the Chen property is trivially equiv- 
alent to the property T{s,t) = T{s,u) (8) T{u,t), implying that T{s,t) = 
r(0, s)*^"^ (8) r(0,t) is a multiplicative increment. In the particular case 
when T{t) = tV, V e M.'^ is a straight line, r(0,t) = expt^f^^ViX* 
belongs to expT^M"^. Easy arguments due to Chow show then that t i-^ 
r(0, is a G-valued path, where g = Lie{G) is the free Lie algebra in 
d generators, generated as a vector space by the successive commutators 
[X*i,X*2], ... In rough path theory, one quotients out 

by ©n>Af+i^"IR'^- Then the quotient Lie algebra Sat is the free N-step 
nilpotent Lie algebra in d generators, and Gn = expg^v is a Carnot group. 
When d = N = 2, Q2 — {X,Y,Z := [X,Y]) is isomorphic to the Heisen- 
berg algebra, and the defect of additivity of the Levy area CAr{s,t) = 
dri{ti) J_*^ dT2{t2) — jl dT2{t2) JI^ dTi{ti), measured by the difference 

CAt{s, t) - CAt{s, u) - CAt{u, t) = 

(ri(t) - V^{u)){T2{u) - V2{s)) - {T2{t) - T2{u)){T^{u) - ri(s)), 

(2.11) 

is encoded into the non-commutativity of the product in the Heisenberg 
group, given by (in the exponential coordinates) (xi,yi,zi) • (x2, 2/2,-22) = 
(xi + yi,X2 + y2,zi + Z2 + \{xiy2 - X2yi)). 

Carnot groups are naturally equipped by homogeneous norms coming 
from the sub-Riemannian Carnot-Caratheodory metric induced by horizon- 
tal geodesies, i.e. minimizing curves with tangent vectors in the Euclidean 
space T^M'^. Then an a-Holder rough path over an a-H61der path P is sim- 
ply an a-Holder GTv-valued path (in geometric terms, an a-Holder section 
of the principal bundle M x Gat) which projects onto P. 

One has the following two fundamental results. 

Proposition 2.2 (see Lyons Lyons-Victoir JJ^J, Friz-Victoir 120^ ) 
Let < < a < 1. 
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1. (Existence theorem) There exists a (highly non-unique) a -rough path 
over any a- Holder path. In geometric terms, one may lift an a- Holder 
section of the quotient bundle M x (Gat/ exp 0^=2^"!^'^) T^^'^ ~ 
M X M'^ into an a~ -Holder section of the principal bundle M x Gn. 

2. (Approximation theorem) Every a- Holder rough path over T is the 
limit in a~ -Holder norm of the corresponding stack of natural iterated 
integrals over some smooth approximation family r^,e — t- ofT. 

The approximation theorem is essential in that it reduces differential 
equations driven by a-Holder paths (through a limiting procedure which is 
often very subtle) to ordinary differential equations. Estimates for solutions 
in a deterministic setting are given in full details in the book by P. Friz and 
N. Victoir (see [20], Chap. 10). 

This general approach is however insufficient for many purposes. Draw- 
backs are of two types: 

- the arguments leading to the existence and approximation theorems 
are abstract, the first theorem relying on the axiom of choice (due to the 
arbitrariness of the lift), and the second one on an interpolation by sub- 
Riemannian geodesies which are notoriously complicated objects; 

- in a random setting, this approach produces in principle deterministic, 
pathwise estimates, which moreover do not depend on the choice of rough 
path. Even in combination with probabilistic tools such as the Malliavin 
calculus, despite beautiful achivements in the case a > 1/4 (such as global 
existence of solutions for bounded potentials [20], existence of a density 
|10j . ergodicity [31]),...) generalizing results known in the case of diffusion 
equations, it does not permit - in the case of stochastic differential equations 
driven by fBm for instance - to produce anything really better than a local 
existence theorem for solutions beyond the barrier a = 1/4. 

Let us mention briefly en passant another related approach due to M. 
Gubinelli |28] and called algebraic integration. Without being too precise, it 
states the existence of a class of T- controlled paths - stable under functional 
transformations and under integration along F, and to which solutions of 
differential equations driven by F belong - whose increments are of the form 

N 

zt-Zs = Y, E C{h,---,in)rl{ii,...,in) (2.12) 

n=l 

for some functions C^{ii, . . . ,in), up to a remainder pts such that pts = 
0{\t — sl^^'^ , with e > 0. The right-hand side of (I2.12p - viewed as a function 
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of i - is a linear combination of the components of the rough path T, while 
the remainder is sufficiently regular so that conventional estimates apply. 
This essentially avoids the use of smooth approximations and requires only 
the knowledge of the quantities F"^, n < N. 

2.3 Fourier normal ordering 

In contrast with this geometric approach, the point of view developed by 
the second author is that a rough path over an irregular path T is something 
"essentially arbitrary", and that one should rather look for explicitly con- 
structed rough paths with "good" properties, which allow better estimates 
than the general ones. 

Let us summarize very roughly the results obtained so far in the following 
Proposition: 

Proposition 2.3 (see JM IM lE^) 

1. A rough path is uniquely determined by an algorithm called Fourier 
normal ordering algorithm from its tree data, which are generalized 
Fourier normal ordered skeleton integrals on domains indexed by trees. 
As a consequence, any arbitrary set of tree data produces a formal 
rough path (see Definition \2.1\) . 

2. Tree data yielding Holder-continuous rough paths by Fourier normal 
ordering may be obtained by various, explicit regularization schemes 
applied to Fourier normal ordered tree skeleton integrals, using multi- 
scale methods and inspired by the renormalization of Feynman graphs. 
In particular, one may construct rough paths B = (B^^, . . . , B^) over 
fBm such that B^^ is in the j-th chaos of fBm 

Fourier normal ordering consists as in section 2 in (1) cutting iterated 
integrals like 1*^(1, . . . ,n) := Jl dTi{ti) £^ dViiti) . . . dr„(t„) into n! 
pieces by applying the Fourier projection operators := ^(l|5^(i)|<...<|.^<^(„)|)i 
where a ranges in the group of permutations of {1, . . . , n}; (2) rewriting each 
piece V^IyO-i ■ ■ ■ ,n) as a Fourier normal ordered integral over the inverse 
image of the simplex {t > ti >...> t„ > s} by cr by using Fubini's theo- 
rem. The inverse image of the simplex decomposes as a union of elementary 

^''i.e. may be written as a j-linear integral expression in terms of B. 
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domains indexed by trees I I. 

Thus the role of Fourier normal ordering is twofold: (1) it allows a 
general algebraic (combinatorial) classification of (formal) rough paths; (2) 
it induces a correct addition of the Holder regularity indices of the tree data 
when recombining them by the Fourier normal ordering algorithm. We have 
seen an example of this when we estimated the variance of the boundary 
terms in section 1. 

The rough paths described in the above Proposition, in the case of fBm, 
say, are not obtained by an explicit limiting procedure; yet they suggest very 
strongly that the construction of rough paths is closely related to renormal- 
ization in quantum field theory. The purpose of the present series of articles 
is to give a probabilistic construction coming directly from quantum field 
theory. We actually conjecture that (some of) the rough paths of the above 
Proposition may be obtained by some limiting procedure from the construc- 
tion of the next sections. 



3 Definition of the interaction 

We recall that /J dBi{ti) f^^ ^-62(^2) represents the area between the straight 
line connecting {Bi{s), B2{s)) to {Bi(t), B2{t)) and the curve. If the curve 
turns right, resp. left, then the Levy area increases, resp. decreases. We 
have seen that represents in some sense the singular part of the Levy 
area. 

It is conceivable that Bi,B2 or 0i,02 represent the idealized, strongly 
self-correlated motion in of a particle, which - although rotation-invariant 
- may not (probably as a consequence of a mechanical or electromagnetic 
rigidity due to the macroscopic dimension of the particle, or any other sim- 
ilar phenomenon) turn absolutely freely. A natural quantum field theoretic 
description of this rigidity phenomenon is to add an interaction Lagrangian 
of the form Cint = (dA^)"^. The fundamental intuition here is that the field 
B is in some sense a mesoscopic field, while A^ depends on microscopic 
details of the theory. 

This is explained in great accuracy in |40j . in a mathematical language. 
A. Lejay shows how a path F may be modified by inserting microscopic 

Given a rooted tree with n vertices indexed by 1, . . . , n, one integrates over the domain 
with coordinates t\,...,tn £ [s, such that ti < tj whenever the vertex i is above the 
vertex j. When the tree is simply a trunk tree with no branching, one gets a usual iterated 
integral of order n. 
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bubbles all along, resulting in the limit in a path which is indistinguish- 
able from the original one, while the Levy area has been corrected by 
an arbitrary amount. Let us give a very simple example. Take for T a 

straight line T{t) ~ q ) ' insert (somewhat artificially) microscopic 

bubbles of size e = M~^p (covered in a time 0{M~p)) at times which 
are multiples of M~p. Then the resulting path r*" has a Levy area of or- 
der MP ■ {M-'^Pf 

^p— i-oo c») while — > r in a -Holder norm whenever 
a- <a since l(r-W-r-(s))-(r(t)-r(.))| ^ ^ 

\t sj" 

verse process of removing microscopic bubbles of a given path so as to make 
its Levy area finite is of course much more hazardous, and looks a little bit 
like an "inverse Joule expansion" (i.e. like putting back all the molecules of 
a gas into the left compartment of a container after removing the wall which 
separated it from the right compartment, a statistical physicist's nightmare, 
sometimes called "Maxwell's devil"). 

Summarizing the above discussion, one must search for an interaction 
which cures the ultra-violet divergences of the microscopic scale, without 
modifying the theory at mesoscopic scale. This is where quantum field the- 
ory comes into play. The interested reader may refer to several excellent 
treatises on the subject (see e.g. [53] or [38] )• It is impossible to give 
here a self-contained introduction to this theory which is one of the main 
foundations of the modern physics of both high-energy particles and con- 
densed matter. Let us however explain in an informal way the most essential 
concepts, and introduce some useful terminology, in order to fill in the gap 
between probability theory and physics. We have tried to make the next two 
definitions as precise and as general as possible. In our case the space-time 
dimension D is simply one. 

Definition 3.1 (ultra-violet cut-off) 1. Let M > I be a constant, 
and x° : ^ resp. 0, non-negative, compactly supported 
function such that x*^ = 1 in a neighbourhood of 0, resp. = {) in a 
neighbourhood of and = 1 in a neighbourhood of the hypersquare 
supj=x,...._D ICjl = 1- These two functions may be chosen such that 
(X"')j>i); '^'i'th x' '■= X^i^^"'^')} define a partition of unity, i.e. 
+ J2j>i X"' = 1- Let p € N. Then the ultra-violet cut-off at scale p 

of a function f-.R^^R'^ts f^P := J-^ ^ [Ej=o^'(0] -^/(O). 
where J- is the Fourier transformation. Roughly speaking, the ultra- 
violet cut-off cuts away Fourier components o/ momentum ^ such that 
1^1 > MP. 
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2. Let C^{x,y) := C<^(x — y) be the covariance of a stationary Gaus- 
sian field (j) '• ~^ Then (f) has same law as the series of 

independent Gaussian fields ^j>o where (p^ has covariance ker- 
nel C-jp := 1-^ . The ultra-violet cut-off at scale 
p of the Gaussian field (p is then (ff^f := Yl]tfo^> '^iih covariance 

Note that (at least for a good choice of the functions x*^ ) ) the Fourier 
transform of ^ is supported on the union of two dyadic slices, M^~^ < 
1^1 < M^+^. In principle one may extend this decomposition to negative 
scale indices j, so that the limit j —> —oo describes the correlations at large 
distances. In our model however - and this makes it very different with 
respect to classical models in statistical physics, see comments below - it 
is only the transition from the microscopic scale p to the mesoscopic scale 
which is non-trivial, and one may essentially restrict to positive indices j. 

Definition 3.2 (interacting fields) Let ip : be a vector-valued 

Gaussian process on M.^ , X := {Xi, . . . , Xq) a set of real parameters, and 
Pi, . . . , Pq (q > 1) homogeneous polynomials on R'^ x (M'^)^. Then the inter- 
acting theory with interaction Lagrangian Cint{4>){x) = ^p=i XpPp{(p{x); V(;/!>(x)) 
is (provided it exists!) the weak limit F(d(p) of the penalized measures 

n,v,pid(f>) := J-e-i'vA„*(<^-''){a.)<ia.rf^^P(^|^)^ (3.1) 

when the volume \V\ and the ultra-violet scale p go to infinity, where: V C 
is compact; dp^P{(j)\y) is the Gaussian measure corresponding to the 
cut-off field 4>^P restricted to the finite volume V; Zv,p is a normalization 
constant called partition function by reference to Gibbs measures. 

In general, (p is stationary, which accounts for the finite volume cut-off 
V, and Jy Cint{(p^^){x)dx diverges when p — >■ oo, which accounts for the 
ultra-violet cut-off at scale p. The parameters Ai,...,Aq are called bare 
coupling constants. Usually the inverse of the covariance kernel of (/> is a 
differential operator of the form C^^ = AyV^ -|- m^, where m is called the 
mass. (In the case of our model, C^^ contains a fractional derivative opera- 
tor instead, but the present discussion remains valid). Formally (forgetting 
about the cut-offs) dp{(f)) gives the trajectories a weight proportional to the 
Onsager-Machlup functional e~2^^^'^'^ +™ )4>,4') ^ gQ ^^j^g parameters Ay and 

play a role similar to the coupling constants Ai, . . . , Ag, and the sum of 
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the interaction Lagrangian and of the Onsager-Machlup functional is cahed 
simply the Lagrangian. 

In general also, (p is self-similar (or at least asymptotically self-similar at 
short distances), so the term in the Lagrangian V(/)(x))(ix has a cer- 

tain degree of homogeneity with respect to a change of scale x ^ ax or equiv- 
alently i— )• a~^^ after a Fourier transform, which gives the main behaviour 
at large momenta ^ - or equivalently at short distances ~ of the correla- 
tions (or so-called n-point correlation functions) ((/)j^(3;i) . . . (/)j„(x„))i/,p := 
/0i^(xi) . (x„)PA,y,p(#). 

Here we take a high-energy physics point of view. Then the bare scale is p; 
in other words, one uses a cut-off at short distances of order -^p^^ 0, 
and wants to understand the behaviour of the correlations at macroscopic 
distances • In principle, the theory is hopelessly divergent in the limit p — )• 
oo if this degree of homogeneity is negative (the so-called non-renormalizable 
case). On the contrary, expanding the exponential e~ ^''^'^^^^''^'1'^^^^'^^ into 
a series leads to only a finite number of diverging terms (called diverging 
Feynman diagrams) if the degree of homogeneity is positive (the so-called 
super-renormalizable case). When this degree of homogeneity is zero (the 
so-called just renormalizable case, often the most interesting one in practice) 
closer inspection is needed. In all cases, for a large variety of models, one 
obtains by iterated integration with respect to highest Fourier scales (i.e. 
with respect to the field components </>^~^, . . . , 0-''''^) an effective theory 
at scale j which may be described in terms of the same Lagrangian but 
with so-called renormalized parameters, by opposition to the bare parame- 
ters, Xp ~^ Xp or Aa Ay, {rn^y . One obtains in general a flow for 
the parameters, i.e. equations of the type (Ay, (m^)-^; {Xp),p' = 1, . . . ,q) := 
F{X^^ , {m?y~^^; {Xp~^),p' = 1, . . . ,q). Solving this flow down to small val- 
ues of J is then the main task of renormalization. An interesting case is when 
one may show that the contribution of the renormalized vertex ApPp((/>, V<j)) 
goes to zero at distances which are large with respect to the bare scale; then 
the theory is said to be asymptotically free at large distances. The best- 
known examples of this behaviour are maybe the weakly avoiding path or 
the (/)^-theory, both in dimension D = A ; see [46^ [22| [TT] for rigorous results 
0. Our model is original for it combines in some sense features of models 

^®In statistical physics, the size of the lattice usually gives an explicit cut-ofT, so one 
may take p = 0. 

^^On the other hand, in high-energy physics, the main example in this respect is that 
of asymptotic freedom at short dtstances (or equivalently high energy) of quarks [301 154) . 
so exactly the opposite point of view with respect to the one we adopt here. 
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of both high-energy physics and statistical physics: namely, the bare scale 
is 0{M~^P), but the theory is asymptotically free at large distances. Letting 
p — >■ oo, the interaction disappears at all finite scales, hence one retrieves in 
the end a Gaussian theory, in which, however, the singular part of the Levy 
area has been cancelled. 

Perturbative methods are by far the most common in physics, because 
they are accessible to non-experts. They rely on an asymptotic analysis 
of the quantities obtained by expanding into power series in the coupling 
constants the exponential weight q- I ^intW{x)dx ^ These are conventionally 
represented as Feynman graphs (we shall show some of these later on for 
our model). Unfortunately, in all interesting cases, the series diverges by 
and large because of huge combinatorial factors, hence perturbative theory 
has only a heuristic status. Constructive methods, on the other hand (when 
they work!), are based on particularly clever finite Taylor expansions, scale 
after scale, and produce converging series (but not power series!); in other 
terms, they are rigorous. However, the technical apparatus needed to explain 
constructive field theory is much more sophisticated. 

Let us now return to the discussion of our model after this long paren- 
thesis. In order to keep track of the degree of homogeneity of the fields - 
and to obtain eventually the expected Holder regularity indices for iterated 
integrals - we need here a just renormalizahle theory (or, in other terms, an 
integrated interaction which is homogeneous of degree 0). Since {dA^Y is 
homogeneous of degree (4a — 2) in time, one shall use in fact a non-local 
interaction lagrangian, \c!^ j j \ti — t2\~'^"Cint{4'i:4'2){ti,t2)dtidt2, where 

Cint (01 , 02) (ti , t2) = { dA+ {ti)dA+ {t2) + dA- {ti)dA- (t2) } , (3.2) 

which is positive for a < 1/4 since the kernel |ti — 12|~^" is locally integrable 
and positive definite. Thus the Gaussian measure is penalized by the sin- 
gular exponential weight e"^ / /^-*('?^i''^2)(ii,t2)|ti-t2r*"diidt2 ^ Equivalently, 
using the so-called Hubbard-Stratonovich transformation 0, we introduce 
two independent exchange particle fields a± = cr±{t) with covariance kernel 
Ca^{s,t) = Ccj_^{t — s) = E(T±(s)(T±(t) = c'J\s — t\~'^°' and rewrite (let- 
ting dp,{4>), resp. dfi{a) be the Gaussian measure associated to (p, resp. 
a = (o"+,cr_)) the partition function Z := Z{X), 

Z := /'e"^-^''":^ l*^"*^l~''"'^™*^'^i''^^^^*^'*^^'^*^'^*^d;u(0) (3.3) 



which is nothing else but an infinite-dimensional extension of the Fourier transform 
E[e'^^] = e-"^^^^"^ for a random variable X ~ J\f{0, a^) 
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as 

Z:= y"e--t^™*(<^i''^2''^)W'^*(i/i(0)d^((T), (3.4) 

where 

C.int{<\>xA2.(y){t) = iA {dA+{t)a+{t) - dA~{t)a.{t)) . (3.5) 

All of this is ill-defined mathematically since (1) cj is a distribution- valued 
process and dA^ is not defined at a// when a < 1/4; (2) one integrates over M 
a translation- invariant quantity (note that (j)i,(j)2,cr are all stationary fields) . 

4 Heuristic perturbative proof of convergence 

Let us now explain the basics of perturbative quantum field theory, and show 
how it suggests (at least heuristically) the assertions of Theorem 0.1. The 

general idea is to expand formally the exponential of the Lagrangian in order 
to compute polynomial moments, ^pl{xl) . . . 'i/;„(x„)e~ -^^'"^^'^i''^^''^^^*)'^* 
also called n-point functions and denoted by ('0i(xi) . . . ^/'„(xn))A! V'j ■ 
01,02,0-+ or as ^ Y.n>Q ^-^^ ■ ■ ■ ipn{xn) (/ £^int{ ■ ;t)cit)"']. We 

do not bother too much about the volume and ultra-violet cut-off here, and 
write ( • )a instead of ( • )a,v,p- Recall first the following classical combina- 
torial facts. A good reference for perturbative expansions in quantum field 
theory is e.g. [38j. 



Proposition 4.1 1. (Wick's formula) Let X = {Xi, . . . ,X2n) be o, (cen- 
tered) Gaussian vector. Then 

E[X,...X2n]= Yl E[Xi,X,,]...E[X,,„_,X,,J, (4.1) 

where the indices range over all pairings of the indices 1, . . . ,2n. Each 
term in the sum is represented as a graph with 2n points connected two 
by two. 

2. (connected moments) Let { ■ ) := E[e^(^)] ^ penalized measure, 
where X = {Xi,X2, ■ ■ ■) is a (centered) Gaussian vector, and ^{X) is a 
polynomial in Xi,X2, .... Then the connected expectation {Xi . . . Xn)c 
(c for connected) is (formally at least) the sum of all connected graphs 
obtained by (i) expanding the exponential; (ii) applying Wick's formula 
and drawing links between the paired points; ( Hi ) identifying all points 
coming from the same vertex, i.e. from the same monomial in ^{X) 
descended from the exponential. 
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Connected expectations exclude in particular vacuum contributions, i.e. 
terms of the form E[e*(^)]E[Xi . . . X„] = ZE[Xi . . . Discarding these 
contributions can be shown to provide automatically the normalizing factor 
^. Then usual expectations {Xi . . . Xn) are obtained by taking all possible 
splittings of {1, ... , n} into disjoint subsets Iitt). . .l±)/p and summing over the 
products of connected expectations ^p'}2h,...jp(Wieh "^^^^^ • • • (DiG/p ^ijc- 
In practice the last operation is trivial for two-point functions {Xi^Xi^) if 
by parity (which is often the case in quantum field theory) the one-point 
functions (Xj)c vanish, so that (XjjXjj) = {Xi-^Xi^)c- 

Let us return to our model. Using a straightforward extension of the 
above Proposition, one may represent {il)i{xi) . . ■ipn{xn))x, tp = 4> or a as 
a sum over Feynman diagrams, ^]^A{T), where F ranges over a set of 
diagrams with n external legs, and ^(F) € M is the evaluation of the cor- 
responding diagram (see examples below); connected expectations will then 
be obtained as a sum over connected Feynman diagrams. More precisely, 
one obtains formally a (diverging) power series in A, Yln>o Ylr„ ^(^n), 
where F„ ranges over the set of Feynman diagrams with n vertices. The 
Gaussian integration by parts formulaic yields a so-called Schwinger- Dyson 
identity, 



{dA^{x)dA^{y))x 



with Fourier transform 



1 



A2Z(A) 
1 



5a+{y) 5a+{x) 



■ J Ci„t{4>l,(t>2,cr+){t)dt 



A2Z(A) "+^''^'(5a+(x) 
— [-C;l{x,y) + ((C->+)(x)(C->+)(y))A 



(4.2) 



^iei'-'"[i-ier-'"(i(-FcT+)(op)A]. (4.3) 



By parity, {\J'{dA^){^)\'^)\ is a power series in A^. 

Introduce an ultra-violet cut-off at scale p as in Definition 13.11 For 
the simplicity of the exposition we shall actually use a brute-force ultra- 
violet cut-off at momentum ilf, i.e. cut off all Fourier components with 



^®an infinite-dimensional extension of the well-known formula for Gaussian vectors, 
E [dx,F(Xi,. . .,X„)] = Ej C''^{i,j)E [XjF{Xi, . . .,X„)] if C is the covariance matrix of 
{Xi,...,X„). 



25 




Figure 1: Bubble diagram witti 2 vertices. By momentum conservation ^ = ^1+^2, 
which leaves out one free internal momentum. 




Figure 2: More complicated bubble diagram with 4 vertices. By momentum con- 
servation = Ci + = + C2 and ^1 = ^' + ^2, which leaves out two independent 
internal momenta. 

momentum |^j > M^. After Fourier transformation, J Cint{-',t)dt becomes 

iA jj^^l^l^ I dCid6'^C5o(6+6+?)->^cr+(^)J'(a</>i)(^i)J'</>2(6), minus a similar 
term involving The square of this expression contributes the following 

)a 

rMP 



term of order 0(A2) to (|7'o-+(^)|2); 



I- MP 

(4.4) 

This is the evaluation of the Feynman diagram represented in Fig. [H ac- 
cording to the following rules. 

Definition 4.2 (Feynman rules) A Feynman diagram in our theory is 
made up of (1) bold lines of type i = 1,2, with momenta • • • evaluated 
as E|J^(^j(^j)P = jriT+^j (^) plain lines of type ±, with momenta ■ ■ 
evaluated as E|J^(T±(^)P = (3) vertices where two plain lines - one 

of each type - and a hold line meet, with a momentum conservation rule, 
^ = ib^i it ^2 (depending on the orientation of the lines). The definition of 
the interaction implies the presence of a further derivation - represented by 
the symbol d on the Feynman diagram - on the resp. 4>2-fidd, and a 
momentum scale restriction < |^2|; resp. > \i2\, o-t vertices involving 
a CJ+-, resp. a --field. The derivation translates into a multiplication by 
resp. i^2 when evaluating the diagram. 
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Figure 3: First tliree terms of the bubble series. The renormalized covariance of 
the cr-field is equal to the sum of the series. 



It is sometimes useful to consider the evaluation of the corresponding am- 
putated Feynman diagram, from which the contribution of the external legs 
has been removed. Here for instance, the evaluation of the amputated Feyn- 
man diagram associated to Fig. [T] is (|^|^~^")^ times the previous expres- 
sions, hence is equivalent to the ^-independent expression — i('A^(M'')^~^° 
when p — 7- oo. It is a diverging negative quantity. (Using the Fourier trun- 
cation of Definition 13.11 onlv changes the constant K.) However, resumming 
formally the bubble series as in Fig. [3] yields, starting from the right-hand 
side of eq. (0]1 . 



n>0 



1 



l-4a 



p\l—ia' 



_|_|^|l-4a 



K\\MPm) 



l-4a 



1 -FirA2(MP/|C|)l-4a 



A2 



1^1 



l-4« 



(4.5) 



On the other hand (see Fig. [3|), the bare cr-covariance |^|i_4a 
replaced with the renormalized covariance 



has been 



1 



1 



1 



l-Aa 



l + K\'^{MP/\i\ 



A-Aa 



ll-4a ' 



(4.6) 



which vanishes in the limit p oo. The essential reason for this is of course 
that the oscillating signs (—1)" in the bubble series evaluation - due to the 
fact that the interaction Lagrangian Cint{<t>i-,(t>2-,o') is purely imaginary - 
result by summing in a huge, virtually infinite denominator. Taking into 
account the possible insertion of (T_-lines between (T+-lines amounts to a 
simple change of the constant K. In physical terms, the interaction in 
lg|iL4Q has been screened by a huge mass term K\'^Mf^^~^°'^ -^p^oo +oo 
(see section 3 for the definition of the mass). More complicated diagrams 
contributing to {\{T(Jj^){£^)\^)\, and involving internal cj-lines as in Fig. [2] 
also vanish when p — > oo. Thus there remains simply: 
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(i-^^^(OPi)A = ^ier'-'". (4.7) 

Hence E\A^{t) - ^^(s)|2 < - s\'^°', as in eq. (fTTSD . 

As for the mixed term {dA^{x)dA'^{y))\, its Fourier transform is given 

by ^1^1^"'^" — i_|_x"A^(A/|g|)i-"t" ' '^^^re K" < K due to the constraints 
on the scales for bubbles of mixed type with one 0"+- and one (T_-leg, which 
vanishes in the limit p oo (note the disappearance of the factor 1 compared 
to eq. ()4.5p . due to the fact that E,a~^ {x)a~ {y) = 0). Thus the covariance 
of the two-component u-field has been renormalized to |g|i-4a\(j^^p ? where 
is a two-by-two positive "mass" matrix with eigenvalues ~ X'^MP^^~^°'\ 

Using eq. (jl.lip . one obtains: 

{27rc^f{A{s,tf)x = {\A+{t)-A+{sf)x + {\A-{t)-A-{s)\^)x 



+K\A-^{s,t)- Ag{s,t)\^ 
cos(t - s)Ol^r'"'"a!e + IE \A-^q{s, t) - Aq{s, t)\ 



4 
A2 

(^iri+i^2)|t-s|'" (4.8) 



for some constants Ki, K2. 

Let us now consider briefly other correlations. For a general discussion 
we need the following easy power-counting lemma: 

Lemma 4.3 (power-counting rules) Let T be a Feyman diagram with 
No- external a-lines, N(f, external (p-lines, and NQtf, external d(j)-lines. Then 
the overall degree of homogeneity (in powers of ^) of the evaluation of the 
corresponding amputated diagram - also called: overall degree of divergence 
- is 1 - 2aN„ + aN^, + (a - 1)Nq^. 

Proof. Let: /o-, resp. I^f,, be the number of internal lines of type o", 
resp. (j) or dcj)] I = /o- + be the total number of internal lines; and 
L = / — y + lbe the number of loops, equal to the number of independent 
momenta (one per internal line, minus one per vertex due to momentum 
conservation, plus one due to overall momentum conservation). Since one a- 
and two 0- lines meet at each vertex, one also has the relations 2Iu + Nu = V, 
and 21^ + + Ng^ = 2V. Now the amputated diagram is homogeneous 
to |^|-(i-4a)^<T-(i+2a)7'^+L+y-Afa^ (counting one derivative per vertex, and 
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minus one derivative per external d(p-leg which is not taken into account in 
the evaluation). Putting all these relations together yields the result. □ 
If a diagram is overall divergent, i.e. if its overall degree of divergence 
is positive, then the diagram diverges (except if by chance the coefficient 
of the term of highest degree in ^ vanishes). On the other hand, the fact 
that a diagram is overall convergent (i.e. its overall degree of divergence is 
negative) does not imply that it is convergent, since it may contain overall 
divergent s?x6-diagrams. One must hence study the behaviour of all possible 
diagrams, with arbitrary external leg structure. 



Figure 4: Higher connected moments of the Levy area. 

The above simple power-counting argument shows that the overall degree 
of divergence of a connected diagram with 2n external a-legs is 1 — 4na. For 
n > 2, this is < 1 — 8a < since a > | by hypothesis, so such diagrams 
are overall convergent. By the above arguments, there remain only the 
connected diagrams in the limit A ^ oo, see Fig. HI whose evaluation is 
independent of A. 

General considerations following from the multi-scale expansions (one 
may refer to [65] for a good, accessible presentation, or to [60] for an ap- 
plication to the Gaussian renormalization of iterated integrals evoked in 
subsection 2.3) show that it is enough to consider the behaviour of diagrams 
whose internal legs have higher (or even: much higher) momentum scale 
than external legs, the so-called dangerous diagrams. Then the momentum 
scale constraint on the vertices coming from Fourier normal ordering implies 
that the external legs of dangerous diagrams may be either of type a or of 
type d<j), but not of type <j). Consider now any diagram whose external struc- 
ture contains external 90-legs. By parity it has at least two such external 
legs, and the previous power-counting rules show that such a diagram is 
always overall convergent. 

Finally, the law of the field (j) is left unchanged by the interaction. 
Namely, all non-trivial diagrams contributing e.g. to {(j)i{x)(j)2{x)) \ involve 
internal u-lines which (as previously "shown") vanish in the limit p oo. 
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On the whole, this is the content of Theorem 0.1. 

The art of constructive field theory is to make the previous speculations 
rigorous. It relies on the following considerations, corresponding to the weak 
points (not to say flaws!) in the above arguments: 

1. While going from eq. (j4.4p to (j4.5p . we have replaced the amputated 
bubble diagram evaluation by its asymptotics when p oo, namely, 
~K\'^MP^^~^'^\ which is simply equal to its evaluation at zero exter- 
nal momentum ^, also called local part. Thus we have actually not 
resummed the whole bubble series, but only the corresponding local 
parts, and observed that this was equivalent to adding a mass term of 
the form i^'A^M^^-*^"^") J \a{x)\^dx to the Lagrangian. 

2. The bubble series is really a terribly diverging geometric series. Renor- 
malization must actually be performed scale by scale. Considering 
only bubble diagrams with momentum in the dyadic slice M^~^ < 
1^1 < leads on the other hand to a converging geometric series 
for A small enough since the term between parentheses in eq. (j4.5p . 

-j^j , is then < 1. This is equivalent to integrating out 

the highest field components {a^,(j)P), as explained in section 3. One 
obtains thus a running mass coefficient m'' of order A^M^^"^")^. The 
procedure must then be iterated by going down the scales step by step. 
Since renormalization reduces the covariance of the u-field, the bound 
on A ensuring convergence does not become worse and worse after each 
step. 

3. We neglected more complicated bubble diagrams as in Fig. [2l Al- 
though these have the same order as the simple bubble diagram of 
Fig. [H as follows from the above power-counting rules, taking into 
consideration all possible bubble diagrams lead to a terribly diverging 
power series in A due to the rapidly increasing number of such dia- 
grams in terms of the number of vertices, with a coefficient roughly of 
order n! in front of A". This divergence is actually due to the accu- 
mulation of vertices in a small region of space of size 0{M~^), where 
j is the momentum scale under consideration. Multi-scale cluster ex- 
pansions in constructive field theory, by considering only partial series 
expansions, avoid this dangerous accumulation process. 

4. By splitting each vertex J C^^{-;x)dx into its different scales, there 
may appear fields c/^g^ , a-' with different scales ji ^ j2 7^ j. Tak- 
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ing this into account in a coherent way in the previous partial series 
expansions lead to complicated combinatorial expressions encoded by 
so-called polymers, which are the main object in use in constructive 
field theory. 

5. In the previous vertex splitting, the field with lowest momentum scale 
(ji) J2 or j, depending on the case) is called low-momentum field. 
Even though the cluster expansion in each momentum scale prevents 
an accumulation of vertices in the same region of space, the compound 
effect of all cluster expansions at all scales produces unavoidably ac- 
cumulations of fields with very low momentum in very large regions of 
space, which is a dangerous problem called domination problem. This 
accounts for the addition of the extra boundary term Cf^^^^.y in the in- 
teraction Lagrangian. Writing out this term and explaining its precise 
form would however take us too far away. 
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